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ABSTRACT 

We used the CSO 10.4 meter telescope to image the 350 /im and 450/zm 
continuum and CO J = 6 — 5 line emission of the IRAS 20126+4104 clump. 
The continuum and line observations show that the clump is isolated over a 4 pc 
region and has a radius of ~ 0.5 pc. Our analysis shows that the clump has a 
radial density profile oc r~^'^ for r < 0.1 pc and has oc r^^'^ for r > 0.1 pc which 
suggests the inner region is infalling, while the infall wave has not yet reached 
the outer region. Assuming temperature gradient of r~° '^^, the power law indices 
become oc r~°'^ for r < 0.1 pc and oc r"^" for r > 0.1 pc. Based on a map of the 
flux ratio of 350/im/450/im, we identify three distinct regions: a bipolar feature 
that coincides with the large scale CO bipolar outflow; a cocoon-like region that 
encases the bipolar feature and has a warm surface; and a cold layer outside of 
the cocoon region. The complex patterns of the flux ratio map indicates that the 
clump is no longer uniform in terms of temperature as well as dust properties. 
The CO emission near the systemic velocity traces the dense clump and the 
outer layer of the clump shows narrow line widths (< 3 km s^^). The clump has 
a velocity gradient of ~ 2 km s~^pc~^, which we interpret as due to rotation of 
the clump, as the equilibrium mass (~ 200 Mq) is comparable to the LTE mass 
obtained from the CO line. Over a scale of ~ 1 pc, the clump rotates in the 
opposite sense with respect to the < 0.03 pc disk associated with the (proto)star. 
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This is one of four objects in high-mass and low-mass star forming regions for 
which a discrepancy between the rotation sense of the envelope and the core has 
been found, suggesting that such a complex kinematics may not be unusual in 
star forming regions. 

Subject headings: ISM: clouds — ISM: jets and outflows — ISM: kinematics and 
dynamics — ISM: lines and bands — stars: formation 



1. INTRODUCTION 

High mass star formation, especially in its earliest phase, is poorly understood because 
of the nature of the objects in this phase. First of all, young high mass stars are deeply 
embedded in their dense natal cloud cores and form and evolve on a short time scale, much 
shorter than that of low mass star formation. Moreover, once a massive star has reached 
the zero age main sequence, its strong UV radiation pressure heavily affects the surrounding 
material, thus making it difficult to trace back the initial conditions of the molecular cloud. 
Finally, high mass stars form in clusters at relatively large distances from us (typically a 
few kpc): it is hence not easy to resolve individual young stellar objects, even for the closest 
massive star forming region, Orion (e.g., Plambeck et al. 1982, Greenhill et al. 1998, Beuther 
et al. 2005), 

The early B type massive (proto)star, IRAS 20126-f4104, is a well-studied object which 
provides us the opportunity to study the early phase of massive star formation in a relatively 
simple configuration. This consists of a luminous IRAS source embedded in a dense core 
and associated with both a bipolar outfiow and a circumstellar Keplerian disk (imaged in 
NH3, CH3CN and C^^S; Zhang et al. 1998, Cesaroni et al. 1997, 1999, 2005). The Cahech 
Millimeter Array (OVRO-MMA) was used to image the large scale (~ 0.5 pc) north-south 



outfiow in detail using the CO J = 1 — transitions (jShepherd et al.ll2000l ). It appears that 
this N-S direction is not preserved on a smaller scale (0.1 pc) where a jet oriented SE-NW 
has been imaged in a variety of tracers, such as shocked H2 (Cesaroni et al. 1997, Shephered 
et al. 2000), NH3(3,3)/(4,4) emission (Zhang et al. 1999), SiO and HCO+ emission (Cesaroni 
et al. 1997, 1999), and scattered light in the near infrared (NIR) continuum (e.g., Ayala et 
al. 1998, Sridharan et al. 2005). The position angle [PA) of the jet is ~ 150°, significantly 
different from that of the large scale molecular outfiow [PA ~ —25°; Shepherd et al. 2000, 
Lebron et al. 2006). Such a discrepancy has been interpreted as precession of the jet/outfiow 
axis (Shephered et al. 2000, Cesaroni et al. 2005). 



While many studies have been done on this object, only few of these have focused on 



- 3 - 



the dense pc-scale clump hosting the (proto)star. It is thus of interest to investigate the 
physical properties and evolution of such a clump and its interaction with the outflow from 
the (proto)star. 

Here we present a study of the dense clump and bipolar outflow based on submillimeter 
continuum and spectroscopic observations taken with the 10.4 meter Leighton telescope at 
the Caltech Submillimeter Observatory (CSO) 0. Note that 1.5 kpc is used as the distance of 
the object to estimate physical parameters in this paper. A discussion regarding the distance 
of this object is given in the Appendix. 



2. OBSERVATIONS AND DATA REDUCTION 

Submillimeter continuum and molecular line observations were carried out towards the 
massive dense clump (Table 1). The dust continuum observations at 350 /im and 450 fim were 
made on 2005 May 10 and 11 using the Submillimeter High Angular Resolution Camera II 
(SHARC II, Dowell et al. 2003; Voellmer et al. 2003) mounted on a Nasmyth focus. Both the 
350yum and 450/im beams have simple Gaussian shapes with half-power beam widths (6'hpbw) 
of 8.7" and 9.8" (see Tabl^. The 384 pixels of the bolometer array cover a region of about 
~ 1' X 2.6' at 350 nm. Pointin g and focus were checked regularly. The pointing drift due 



to temperature (jShinnagal |2004| ) was subtracted in the data reduction process. T he final 



point ing accuracy is within ~ 2". The Dish Surface Optimization System (DSOS) (ILeong 



20051 ) was used during the observations: this allows us to apply a gravitational deformation 
correction to the main mirror physically and is thus a key component for high sensitivity 
short submillimeter wave observations. Zenith opacities T225 ~0.06 at the wavelength of 1300 
Hm were measured with the CSO tau meter during the 350/im observations, while during 
the acquisition of the 450/im data we measured T225 ~0.09. The data were reduced with the 
software Comprehensive Reduction Utility (CRUSH, Kovacs 2006) version 1.41. The flux 
was calibrated with Neptune for both wavelengths, with an estimated calibration uncertainty 
of 7% at 350yum and 5% at 450yum. The estimated flux of the planet was 94.6 Jy at 350/im 
on 2005 May 10 and 55.9 Jy at 450/im on 2005 May 11. Note that the angular size of the 
planet was much smaller than the beam sizes at these wavelengths. An area of ~ 9' centered 
at the source position was mapped with the instrument, in box-scan mode. The RMS noise 
of the final maps is about 200 mJy/beam for both bands. After calibration, data analysis 
was done using the program GRAPHIC of the GILDAS software. 

Spectroscopic observations in the CO J = 6 — 5 line (z/= 691.473076 GHz, Goldsmith et 
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al. 1981) were made on 2004 July 8 and 9 using a 690 GHz heterodyne receiver. An area of 
2/7 X 2/7 was sampled with steps of 5" and 20". while the innermost region of 1' 20" x 1' was 
sampled with steps of 5". r225 was between 0.07 and 0.085 during the observations. The SIS 
receiver operation at 4K produced typical single sideband system temperatures (measured 
with the 1.5 GHz bandwidth Acousto-optic spectrometer) between ~ 7000 K and 10900 K 
at 691 GHz, for elevations in the range of ~ 50-70 degrees. The antenna temperature 
calibration was done by the standard chopper wheel method every 30 minutes to 1 hour. 
The data were converted to the main beam brightness temperature (T^h) scale. The ^hpbw 
at 691 GHz is ~10". The main beam efficiency at the frequency of the CO 6 — 5 transition at 
an elevation of 49° was estimated to be ~ 40 %, through observations of Uranus. Pointing 
was checked approximately every two hours. The pointing accuracy is estimated to be ~ 5". 
Daily gain variations as well as gain curve over different elevation angles were derived and 
corrected using CO 6 — 5 spectrum of NGC 7027. Frequency calibration was done every 30 
minutes. The data were reduced using the program CLASS of the GILDAS software. 



RESULTS AND ANALYSIS 



3.1. Submillimeter Dust Continuum 



3.1.1. Overall Structures and Flux Measurements 



Submillimeter continuum observations over a 4 pc region confirmed that the object is 
an isolated massive dense clump (Figured]). The submillimeter dust peak positions in both 
bands coincide with the 1mm and 3mm continuum peak position (RA = 2 0^ 14™ 26.^03, 
Decl. = 41 °13^32."7 (J2000 ) : 1 = 7 8°0r20."2, b = 03°37'58."3)), measured bv lCesaroni ^TaL 



(]1999[ ) and lShepherd et al.l (l2000f ). 



The 350/im image traces a north-south elongated single dense clump with a size of 
~1.3 X 0.7 pc, which fits the silhouette against diffuse ionized emission. The 450/im dust 
continuum, which has a size of about 0.9 x 0.4 pc, looks more compact, with some extension 
towards the south. The clump has a prolate spheroidal shape. The effective radius of the 
clump is defined as R^s = {A/n — 6'hpbw/4)^^^5 where A is the area of 2a contour traced by 
350 fj,m and 450/im continuum. We find -Refrsso = 0.48 pc and ReM5o = 0.37 pc, corresponding 
to angular radii of 66" and 51", respectively. 

The peak flux values at 350/xm and 450/im are 59. 1± 4.1 Jy beam~^ and 27.4±1.2 Jy 
beam~^, respectively. Within the u ncertainties, these valu es are consistent with previous flux 
measurements done with SCUBA (ICesaroni et al.l Il999l ) with 6'hpbw=15", but the 350/im 
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peak flux value is higher than another SHARC measurement at the same wavelength, 34 
Jy beam~^ (lHunter et al.ll200Cll ). with 6'hpbw= H"- All peak flux densities at different wave- 
lengths reported to date are summarized in Tabled Our total flux densities computed inside 
the 2(7 contour are 477 ± 13 Jy at 350/im and 137 ± 14Jy at 450/im, respectively. 



3.1.2. Radial Density Profile 

We approximate the column density (A^(H2)) as a function of clump radius (r) as a 
power law, A^(H2)oc r^. This reflects the density structure of the clump, which to a flrst 
order approximation should also scale as a power law, with index p — 1. Figure [2] represents 
the azimuthally averaged column density proflle derived from the 350/im data. The density 
proflle is plotted from the clump center, i.e. the peak position at 350/xm. A discussion on 
the effect of a temperature gradient on the estimate of the column density proflle will be 
given in § 14. 1[ 

One can see that the inner part of the dense clump has a flatter index compared to that 
of outer part. The power law for 0.03 ^ t < 0.13 pc (corresponding to 4.4" and 17") is 
j,-o.26±o.o3_ j^qIq tliSit this region is larger than ^hpbw35o/2 = 4. "3, so that the flatter index 
observed in the inner region is not due to beam dilution. In the range 0.13 ^ r < 0.3 pc (i.e. 
between 17 and 40"), the column density decreases as Finally, for 0.3 ^ t < 0.7 

pc (40 - 100"), the index p is —1.98 ± 0.02. Steepening of the index in this region is likely 
due to geometrical effects close to the clump border at ~50" from the center. Instead, in the 
two inner regions the estimates obtained are reliable and imply power law indices of —1.26 
and —2.32 for the volume density. These issues will be further discussed in § 14. 1[ 



3.1.3. 350fim/450fim Flux Ratio Map 

The dust continuum maps taken at two different wavelengths allow us to calculate a 
flux ratio map. Note that such a ratio depends on two parameters: the emissivity spectral 
index, /5, and the temperature of the dust grains, T^ust- 

Figure [3] shows a map of the 350/xm/450/im flux ratio. Note that the 350/im data have 
been smoothed to the same angular resolution as the 450/im image. Only data above 2 a 
have been used in the computation of the ratio. The plotted ratio range is between 2 and 3 
in order to show the change of the ratio clearly. 

The flux ratio map reveals a few intriguing features. First of all, the flux ratio towards 
the clump center is low (~ 2) and shows a bipolar structure elongated in the north-south 
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direction. Second, the flux ratio increases towards the outer region and reaches the maximum 
(~ 3) at r ~ 0.3 pc. Such a high flux ratio is found in all directions and traces a region 
elongated in the north-south direction. Third, the flux ratio decreases down to < 2 beyond 
the high flux ratio region at r > 0.4 pc. 

This variation of the flux ratio can be found also in Figure |H Here, the flux ratio 
is calculated over circular annuli with 10" width and 10" step. The error bars are the 
standard deviation of the ratio values measured in each annulus. Note that error of the flux 
ratio increases towards the clump edge, as the intensities of the continuum emission decrease. 

Considering the fact that there is only one prominent heating source, the massive 
(proto)star, inside the dense clump, the flux ratio pattern is expected to be spherically 
symmetry. However, the observed pattern of the flux ratio map looks significantly elon- 
gated, which may be due to the bipolar outflow from the massive (proto)star. This will be 
discussed in detail in § HI 



3.2. Molecular Components Traced with CO Emission 

We use our maps of the clump in the CO J = 6 — 5 line to investigate the physical 
properties of the molecular gas (§ 13.2. ip as well as its low-velocity (§ 13.2.21) and high- velocity 
components (§ 13.2.31) . 



3.2.1. Physical Properties Estimated From CO J = 6 — 5 and J = 7 — 6 Spectra 

In Figure [5] (a), the CO J = 7 — 6 spectrum (Kawamura et al 1999) is compared with 
the CO J = 6 — 5 spectrum to infer the physical properties of the molecular gas. Figure [S](b) 
shows the ratio of the main beam temperatures of t hese transitions. Since the temperature 



scale of the CO 7 — 6 spectrum is a lower limit (IKawamura et al.l 119991 ) , the brightness 
temperature ratio is to be taken as an upper limit. One can see that the red shifted wing 
has a higher ratio than the blue shifted one. 

When dust couples well with molecular gas, the excitation temperature of the gas must 
be close to the temperature estimated from the spectral energy distribution (SED) of the 
dust continuum. Therefore line emission is expected to be optically thin when Tmb is lower 
than Tdust=T'ex (~ 40 K, Shepherd et al. 2000; see also § 14. ip . Thus, the CO emission in the 
line wings is likely optically thin, because T^b^ 40 K at high velocities. 



We now wish to obtain an estimate of the gas temperature in the high velocity gas (i.e. 
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in the line wings) using the result in Figure [5] (b). For this purpose we express the ratio 
between the main beam brightness temperatures of the CO 7 — 6 and 6 — 5 lines as 

2 



-mb CO 7-6 



mb CO 6-5 



/6-5 Ki^r^e) 5^6-5 (^cx) - 5^6-5 (^bg) l-exp(-r6_5) 



where Tmb is the main beam brightness temperature of the CO transition, / is the corre- 
sponding beam filling factor, and Bi,j_^j_^^{T) = ■ [exp( ^'''^~y"^' ) — 1]^^ is the Planck 
function. Here h and k are the Planck and Boltzmann's constants. Tbg is the cosmic mi- 
crowave background temperature (2.7 K). At the frequencies of the CO 7 — 6 and 6 — 5 
transitions {i^7-e= 806.652 GHz; z/6_5= 691.473 GHz), 5^,^ .-(Tbg) and B^^_^iTi,g) are negli- 
gible. As already pointed out, in the line wings r <^ 1, so that [l-exp(-r)] ~ r. Since the 
source is resolved, we also assume ~ 1. Under these approximations, the above equation 
can be re-written as 



Tmh CO 7-6 ^ ^^7-6 ^ y fcT^x J Tj^q 
Tmh CO 6-5 ^^6-5 gxp ( \ - 1 ^6-5 



Under local thermal equilibrium (LTE) condition, the ratio of the optical depths of the two 
transitions can be calculated from the following equation; 

r7-6 . , ^^7-6 ( hu7-6 \ ^^P(Sf) ~ ^ /ON 



r6-5 1^6-5 V ^^ex J exp(^) - 1 

Substituting Equation (3) in the Equation (2), replacing hv-j^^/k = 38.714 K and ui^q/uq^^ = 
7/6, and solving with respect to Tex, one finds 

38.714 K 

= z ^ . (4) 



-ex r o rn 

^ly T^h CO 



In 



Tmb 



Based on this equation, the blue- and red-shifted components have temperatures equal 
to 130 K and 30 K respectively and we thus suggest that the temperature of the blue-shifted 
gas is higher than that of the red-shifted gas. Please note that the difference on the ratio of 
the brightness temperatures of the CO lines could be partly due to difference in densities in 
the two outflow lobes. 

We also attempted an estimate of the optical depth of these transitions using the large 
velocity gradient (LVG) model. Assuming an H2 volume density of 10^—10^ cm~^, a CO 
column density of 10^^ — 10^^ cm~^, an abundance ratio between molecular hydrogen and CO 
of 10^ (see also § I3.2.3p . and a measured line width of 25 km s~^, we obtain tj_q = 1.2 and 
re_5 = 0.7 from the model. The tau ratio is 0.7, similar to that obtainable from Equation (3). 
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3.2.2. Narrow Line-Width Rotating Massive Dense Clump 

Figure [6] shows a comparison between the 350yum dust continuum and CO 6 — 5 hne 
emission with LSR velocity in the range from —6.75 to —0.25 km s~^. This velocity range is 
chosen so that the emission in the range corresponds to the extended components. 

As shown in Figure El the warm gas traced by the CO 6 — 5 line shows good spatial 
correlation with the 350yum dust continuum emission. The molecular component traced by 
CO 6 — 5 is not spherically symmetric. The effective radius of the CO 6 — 5 emission, 
-Reff C06-5; IS 0.56 pc (77"), larger than R^q 350 (§ 13.1.11) . There is a strong CO 6 — 5 peak 
towards the clump center (i.e. the continuum peak). 

The column density of the CO molecule using J = 6 — 5 transition can be derived from 
(e.g.. Equation (Al) of Scoville et al. 1986): 



3fc Te. + ^ jT^^.dv 

87r35co/X^O ■ 6 \_ ^^^^ ' Jc06-5 - Jbg6-5 ' ^ ^ 

where Bqq is the rotational constant of the molecule (57897.5 MHz), z/coe-s the frequency 

'"^C06-5 

of the transition, J(T) = hucoe-^a specific intensity, and v is the velocity. Note that 

Tm.h= i_e-r^ '^b, whcrc Tb is the brightness temperature. For J T^^^dv, we use an integrated 
spectrum of the area inside the 2 sigma contour. Since CO 7—6 data are available only for 
the central region and we don't have isotopic hne data, we assume that the column density 
of molecular hydrogen is proportional to the observed velocity integrated ^^CO emission, 

/ Trahdv. 

The LTE mass of the clump, Mlte, is calculated as follows: 

Mlte = A ■ A^'coe-s ■ ■ Ate ■ "^Ha (6) 

with abundance ratio of molecular hydrogen to CO (10*^), hg mean atomic weight of 
the gas (1.36), and mass the hydrogen molecule. The derived LTE mass is 390 Mq for 
Tcx=40K. Hereafter, 40 K will be used as the temperature of the clump. Please note that 
the optical depth of CO 6-5 line that traces the extended clump components towards limited 
central region is likely to be thick. Accordingly, the LTE mass derived above is lower limit. 

Figure [7] shows the LSR velocities towards selected positions where the full width at 
half maximum (FWHM) of the line is < 3 km s^^. While the east part of the clump has 
LSR velocities ranging between ~ —1 to —3.5 km s~^, in the west side the velocity ranges 
from —3.0 to —5.2 km s~^. The velocity gradient measured over 1 pc is ~ 2 km s~^pc~^. 
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The rotation angular velocity, fi, is ~ 7 x 10~^^ s~^. In Figure [8] we show template spectra 
of the narrow line-width components observed in the east side and west side of the clump. 
This comparison demonstrates the clear velocity shift from east to west observed across 
the clump. Closely looking at the VLSRchange, the velocity shift might be from northeast 
to southwest. The direction of the velocity gradient is almost perpendicular to the major 
axis (directed north-south) of the massive dense clump. Assuming that the velocity shift is 
caused by rotation of the clump and that centrifugal forces balance gravitational ones, one 
can estimate the equilibrium mass of the molecular dense clump as Mj-^t — ^^q '^ = 200 Mq, 
where G is the gravitational constant. This value agrees reasonably well with the LTE mass, 
supporting the idea that the velocity gradient may be largely due to rotation of the clump. 

We note that rotation of the narrow line-width molecular gas component occurs in the 
opposite sense with respect to rotation of the circumstellar Keplerian disk associated with 
the (proto)star. A similar discrepancy in the se nse of rotatio n on different spat ial scales, has 



been observed also in W3 (IHayashi et al.lll989l ) and Orion (IHarris et al.lll983l ). for massive 



star forming regions, and in a starless core, L1521F (MC27) in the Taurus Molecular Cloud 



( IShinnaga et al.ll2004l ). for low-mass star forming regions. This suggests that such a complex 



kinematical structure may not be unusual in the star formation process. 

The measured Avfwhm averaged over the entire CO 6 — 5 emitting region is 3.9 ± 0.2 
km s~^. From the line width, one can derive the virial mass, Mym, of the clump: 



J, J 5 + 2e _Rcfr co6-5 2 



where the volume density is proportional to r~'^. The virial mass is estimated to be 4OOM0for 
this dense clump. One can assess the equilibrium of the dense clump by comparing this value 
with the LTE mass and the one obtained assuming centrifugal equilibrium. We note that 
Mlte — Mj-ot and at the same time, Mlte — ^vm- This suggests that the clump is in an 
equilibrium state and rotation of the core contributes significantly to the equilibrium of the 
clump. 

Figure M shows the non-thermal velocity dispersion, cr^^^, of the narrow line- width 
molecular gas components. Note that the diagram of the velocity width is plotted under 
an assumption of constant temperature over the core. Since the optical depth of the CO 
6 — 5 line is relatively small (see § I3.2.ip . one can estimate the non-thermal velocity dis- 
persion (o"^,^J by subtracting the thermal line width (^^)^^^ from total velocity dispersion, 
(^v,,tai = [<bs - ^^k(^ - ^)]^/^ where a^^,^ is the observed velocity dispersion, m^^^ is 
the average mass unit, and mco = 28 amu is the mass of the CO molecule, o'v^^^^ can be ex- 
pressed as At>FWHM(4 ln2)~^/^, where Awfwhm is the observed hne FWHM. From Figure [TOl 
one can see that the non-thermal velocity dispersion decreases with increasing distance from 
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the clump center. To derive the non-thermal velocity dispersion, 40 K has been used as gas 
kinetic temperature. The dotted line in the diagram shows the thermal velocity dispersion 
when Tk= 40 K. One can see that the thermal velocity dispersion becomes dominant at r= 
0.6 pc (~ 80"). 

3.2.3. CO 6 — 5 Large Scale Bipolar Outflow 

Figure [11] shows the CO J = 6 — 5 spectra taken towards three positions and the mean 
spectrum over a < 120" area around the clump center. One can see prominent wings in 
the spectra taken towards each outflow lobe. Blue- and red-shifted emission is detected over 
the velocity ranges -34.0 < 1/lsr km s-^< -6.0 and 0.0 < VtsR km s-^< +30, respectively. 
The spectrum observed towards the clump center (0", 0") and the mean spectrum clearly 
show both red- and blue-shifted emission. The spectrum taken towards the center position 
has a strong peak at Vlsr=— 3 km s~^ 

Figure [12] shows the distribution of CO J = 6 — 5 blue-shifted and red-shifted emission 
superposed on the 350/im dust continuum emission map. Both outflow lobes are located 
inside the 350/im dust continuum component. The red-shifted bipolar outflow lobe is more 
than twice extended compared with the blue-shifted outflow lobe. Integrated intensities 
of the blue- and red-shifted outflow lobes are measured to be (9.6 ± 0.9) x 10^ and > 
(2.5 ±0.1) X 10^ K km s~^, respectively. The masses of the blue and red shifted outflows are 
estimated to be 28 and 13 Mq, respectively. From the outflow velocity, ~30 km s~^, and the 
extension of the outflow, ~0.6 pc, one can estimate the dynamical time scale of the bipolar 
outflow. The averaged dynamical time scale of the outflow is estimated to be > 1.7 x 10'* 
years. Note that the estimated age is lower limit because the south edge of the red shifted 
outflow lobe isn't fully covered and the inclination angle is unknown. The inclination angle 
affects the size and velocity. The outflow parameters derived from the CO 6 — 5 and 1-0 
lines are summarized in Table [3] The outflow velocities obtained from the 6 — 5 line are 
higher than those from the 1 — 0, while the mean dynamical time scale of the 6 — 5 flow is 
shorter than that of the 1 — 0. A more detailed discussion of the outflow properties will be 
given in § m 
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DISCUSSION 



4.1. Dust Property and Column Density Profile 



To derive the dust mass M^ust, we have used the expression (see Hildebrand 1983), 



M, 



dust 



dust ) 



where F^, is the measured flux density, D the distance to the object, and i?;y(Tdust) the Planck 
function. The coefficient ( ku is dust opacity), is estimated to be 0.1 g- cm"^ at 250/im 
in NGC7023 (IHildebrandl Il983l ). depends on the frequency as n^, cc u^. (3 may range 
from ~ 1 at A < 200/im to > 2 at ~1 mm (Hildebrand 1983 and the references therein). 
Extrapolating from 250/im to 350/im, one obtains i^sbOGHz — 5.6 cm^ g~^. /3 can be 
estimated from the flux ratio of 350/xm and 450/im using the equation 



^2 



(Tc 



dust I 



f2 



dust ) 



(9) 



where VLy is the beam solid angle. Note that in our case ^ ^ 1, because the 350/im image 
has been smoothed to the same angular resolution as the 450yU.m image. 

Assuming that the clump mass equals M^ot — 200 M©, one may use Equation ([8]) 
to derive an estimate of the mean clump temperature, T^ust = 40 K. Then Equation ([9]) 
gives (3 ~ 1.7. This value lies within the typical range 1.5 ^ /3 < 2.4 for high mass star 
forming regions (e.g., Molinari et al. 2000). The dust temperatu re is more or les s cons istent 
with the temperature, ~ 44 K, deri ved from the SEP fitt ing by lShepherd et al.l (120001 ) . but 
lower than the 60 K estimated by ICesaroni et al.l (119991 ). Note that the errors on j3 and 
especi ally T^ust may be s ignificant, as Mrot (200 Mq) might not represent the true clump 



mass. 



Hofner et al.l (120071 ) estimates a clump mass of ~ 400 M0from their SED fitting. The 



column density towards the clump center is estimated to be ~ 3.4 x 10 cm~ , using the 
values of (3 and T^ust described above. This implies a volume density of ~ 10^ cm~'^ towards 
the clump center. Towards the outflow lobes, the estimated number density becomes ~ 10^ 
cm~^. 

A temperature gradient is likely to exist inside the clump, as this is internally heated 
by the embedded (proto)star. Such a gradient is likely to be similar to a power law of the 
type T(r) oc r~'^/^^'^^^ (Doty & Leung 1994), with r distance from the clump center. In 
our case (3 = 1.7 and hence the the power law index is —0.35. The effect of this gradient 
on the estimated density slopes inside the clump is to change the power law indices of the 



inner (0.13 



< 



< 



0.3 pc) and outer (0.3 



< 



< 



0.7 pc) regions to ~ —0.9 and 
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respectively. This has been assessed by recalculating the column density profile assuming 
the clump is spherically symmetric. These indices are not significantly different from those 
derived without taking the temperature gradient into account. For the non-thermal velocity 
dispersion described in § 13. 2. 2^ the velocity dispersion won't change much (order of 0.01 km 
s~^) even when we consider an temperature gradient of T(r)oc r~°'^^. 

An index of —1.26 for the 0.03 ^ r < 0.13 pc region is only slightly shallower than 

3 

that typical of a region undergoing free-fall {p{r) oc r~2) - as in the runaway collapse model 
(Larson-Penston-Hun ter solution, Lars on 1969; Penston 1969; Hunter 1977) and the inside- 



out collapse model by lShu et al.l (119871 ). Instead, an index of —2.32 for the 0.13 < r < 0.3 pc 
region resembles the density profile p(r) oc of a system in hydrostatic equilibrium. We 
thus speculate that the inner region of the clump is experiencing infall, whereas the "infalling 
wave" has not yet reached the outer region beyond r > 0.1 pc. 



4.2. Internal Structures of The Dense Clump 

4-2.1. CO 6 — 5 Outflow and Shocked H2 Knots 

From Figure [T2|, one can see a spatial coincidence between the outfiow lobes and the 
shocked H2 knots. All six knots lie inside the CO 6—5 bipolar outfiow lobes. More specifically, 
the CO 6 — 5 emission shows clear peaks close to the H2-5 and H2-6 knots located in the red 
shifted lobe. Such a correlation between CO 6 — 5 peaks and shocked H2 knots indicates that 
they are physically associated. Note that, unlike the CO 6 — 5 emission, no peak towards 
the knot positions is found in the CO 1-0 line, as one can see in Figure 1 of Shepherd et 
al. (2000). A likely interpretation is that CO 6 — 5 traces more excited gas compared to 
the 1 — transition. This might be due to the gas at the knot locations being heated and 
compressed by interaction with the jet powered by the (proto)star. 



4.2.2. Implications from the Suhmillimeter Dust Radio Map 

There are three major components in the submillimeter fiux ratio map (see § 13.1.31 and 
Figure [3]), namely: (1) a bipolar feature in the central region of the clump, with fiux ratio 
~ 2— 2.2; (2) a region ~ 0.5 pc in diameter and elongated north-south, with fiux ratio 
starting from ~ 2 — 3 and increasing upto > 3 close to the outer border; and (3) region of 
lower fiux ratio, < 2, outside of component (2). Component (1) is located inside component 
(2). The PA of the bipolar structure is about —17°. Note that features (1) and (2) are 
detected above a 3(j level in the 450 yum map and above a 5(j level in the 350 /im map. Note 
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that 450/im emission in component (3) is detected at a 2 — 3cr level, whereas the 350/im 
emission is much stronger. Therefore, we believe that all the three features represent real 
structures. 

Equation (Q shows that the ratio between the 350 and 450/im fluxes is an increasing 
function of both P and (albeit more weakly) Tdust- It is thus impossible to decide whether 
a change in the ratio is due to a variation of (3, temperature, or both. For example, for 
T = 40 K, a /3 of 1 and 2 implies a flux ratio of 1.9 and 2.4, respectively. Vice versa, 
for (3=1.7, a temperature of 20, 40, and 60 K implies a flux ratio of 1.9, 2.3, and 2.4, 
respectively. From these numbers one can see that, in order to reproduce flux ratios as large 
(~3) as those measured on the surface of component (2), very likely an increase of both (3 
and Tdust is needed. Hence component (2), which encases the bipolar feature, might be a 
cocoon of dense material with warm surface. On the other hand, the lower flux ratio value 
of component (3) indicates that Tdust and/or P in the region ahead of the outflow lobes 
are lower compared to component (2). We speculate that the higher Tdust and (3 observed 
on the surface of the cocoon represented by component (2) might be due to a weak shock, 
propagating from the (proto)star and possibly created by the outflow. Considering a sound 
speed of ~ 0.8 km s~^, the crossing time of the weak shock would be < 3 x 10^ years. Note 
that this is likely an upper hmit to the age of the outflow (and perhaps of the (proto)star), 
as shock waves may propagate with larger velocities than the sound speed. 

How to explain the low flux ratio in component (1)? The fact that the shape of this 
region matches very well the outflow lobes suggests a connection between the two. This 
suggests that the flow might have swept away the smallest grains from that region, thus 
biasing the grain size towards the largest ones. This phenomenon would have the effect of 
reducing the slope of the spectral index P and hence of the flux ratio. 

In conclusion, the flux ratio variations across the clump strongly suggest that the mas- 
sive dense clump is no longer uniform in terms of temperature as well as dust properties. 
To construct a model to explain the observed features, one must keep in mind such inhomo- 
geneities. 

4-2.3. Comparison Between Suhmillimeter Flux Ratio Map And Molecular Components 

Here we make a detailed comparison between the flux ratio map and the structure of 
the dense clump and outflow. Figure [13] shows the CO 6 — 5 and CO 1 — outflow lobes 
superposed on the flux ratio map. We have already discussed the good correlation between 
the CO outflow and the bipolar component (1) of the flux ratio map. Now, we note also 



-14- 



that the peak of the blue shifted CO 6-5 lobe is located on the west edge of the northern 
bipolar feature, which may be due to interaction with the jet detected in SiO and H2 line 
emission (Cesaroni et al. 1999) and oriented with PA~ —60°. This jet might have excited 
the molecular gas lying on the west side of the northern lobe. Nothing similar is seen to the 
SW in correspondence to the CO 6 — 5 red shifted outflow lobe. There are two distinct CO 
6 — 5 red shifted peaks associated with the shocked H2 knots, H2-5 and -6, located at (14", 
-30") and (17", -48") respectively, as described in § 14.2. 1[ At the position of the H2-5 knot, 
a high value of the flux ratio (~ 3) is observed, suggesting relatively large values for both 
the temperature and (3. In correspondence to the CO 6 — 5 peak associated with knot H2-6, 
the flux ratio is < 2.5. The difference in flux ratio between H2-5 and H2-6 may indicate that 
the temperature and/or (3 at the position of knot H2-6 could be lower than those of H2-5. 

5. SUMMARY 

An observational study of the IRAS 20126+4104 clump is presented. The main findings 
of this study are summarized as follows. 

(I) The 350/im and 450/im dust continuum images revealed that the power law index of 
the volume density in the inner part of the clump is — 1.26±0.03 (r < 0.1 pc), shallower than 
that in the outer part, —2.32 ± 0.02 (0.3 < r < 0.1 pc). The power law indices suggest that 
the inner region might be infalling, while the outer region of the clump is still in hydrostatic 
equilibrium. 

(II) A map of the peak velocity of narrow line- width CO 6-5 components reveals a veloc- 
ity gradient across the clump, suggesting rotation about an axis oriented northwest— southeast 
The rotation occurs in the opposite sense with respect to the circumstellar Keplerian disk 
associated with the massive (proto)star. 

(III) Assuming centrifugal equilibrium, one estimates an enclosed mass M^^t of ~ 
2OOM0. Estimated Mlte and MyiR (~ 400Mq) of the clump are comparable to M^^t- 
This indicates that the clump is in equilibrium state as a whole and rotation contributes a 
large fraction to support the clump from gravitational collapse on large clump scale. At the 
radius of 0.1 pc, the rotation may be disconnected from the inner region and the medium 
inside of r ~ 0.1 pc may be infalling. This might be partly related to the disagreement 
of the rotational axes of the large scale clump and of the disk associated with the massive 
(proto)star. 

(IV) The mean spectral index (3 and dust temperature are estimated to be 1.7 and 40 K, 
respectively. 
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(V) The 350/im/450/im ratio map outlines three major components, namely: (1) a 
bipolar feature, which may represent the region swept by the bipolar outflow lobes; (2) a 
dense cocoon ~ 0.5 pc in diameter, with a warm surface; and (3) a low temperature/low (3 
component outside component (2). These three components indicate that the temperature as 
well as the nature of the dust (grain size distribution and composition) may vary significantly 
inside the dense clump. The warm layer of component (2) might have been created by a 
shock propagating from the (proto)star. The crossing time of the shock is estimated to be 
> 3 X 10^ years. 

(VI) Strong correlation between CO 6 — 5 peaks and two of the shocked H2 knots are 
identified, indicating they are physically associated. 
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A. Appendix: Distance To The Massive Dense Clump 



IRAS 20126+4104 is located in the no rthwest side of the Cygnus X region. This region 
( iPiddington fc Minnettlll952l : lDavieslll957l ) is one of the most active star forming complexes 
in the Galaxy. Here one can find a large variety of objects, including OB associations, atomic 
diffuse clouds, molecular dense clouds/cores, evolved stars, and supernova remnants. The 
distance of the massive dense clump under study is not well established as it lies on the tan- 
gential direction of the spiral arm or bridg e structure to which the Sun belongs in our Galaxy 
( 1Sharplesslll965l : iKerr fc Westerhoutlll965l ). Many studies reporting on o bservations of IRAS 
20126 +4104 ado pt a distance of 1.7 kpc, b ased on the value quoted by lDame fc Thaddeus 
( 119851 ). However, iDame fc Thaddeud (119851 ) do not explain how such a distance was derived. 
The kinematical estimate using a galactic rotation model may not be accurate for this source 
because of the location of the source relative to the Galactic Center. An alternative method 
is that of looking for a possible interaction of the IRAS 20126+4104 clump with nearby 
objects whose distances are known. From the north-east to the north of the source, there is 
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a chain of Hll regions, whose northern edge coincides with IC1318a, the biggest HII region 
in the cluster. IRAS 20126+4104 is located close to the south-western edge of the chain. 
The angular distance from the closest HII region is about 7'. In addition, there is a num- 
ber of HII regions to the east of the sou rce and one supernova remnant (SNR), G 78.2+2.1 



(IHiggs et al.lll977l : lLandecker et al.lllQSOl ) with diameter of ~62" to the south-east of it. The 
angular distance between the dense clump and the supernova remnant is about 1.5°, along 
a direction with PA ~ 120°. To the south-east of the object, there is an O star, HD193322 
(09V). The angular separation between HD193322 and IRAS 20126+4 W4js_abou^_5]/jnie 



distance of the O star is estimated to be ~1 kpc from photometry (jCruz-Gonzalez et al. 



19741 ). To the south-west, there is no strong HII region. Ionized gas appears to surround the 
massive dense clump except on the south-west side, as shown in Figure 1 (§ 13. ip . which is 
evidence of interaction between our clump and all these objects. The ionized gas is likely due 
to irradiation from these HII regions and from the SNR. As the distance of these sources is 
reported to be 1.5 ± 0.5 kpc (Iksanov 1960, Dickel et al. 1969, Higgs et al. 1977, Landecker 
et al. 1980, Piepenbrink & Wendker 1988, and references therein), in this paper we assume 
a distance of 1.5 kpc for IRAS 20126+4104. 
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Table 1. Observations 



Emission 


Frequency 


Bandwidth 


Telescope / 


©HPBW^ 


FOV*^ 


Velocity resolution 


Reference 




(GHz) 


(MHz) 


Instrument 


(") 


(arcmin^) 


(kms~^) 




350/im continuum 


849.4 


102900 


CSO/SHARC H 


8.7 


9x9 




This work 


450/im continuum 


658.3 


67900 


CSO/SHARC H 


9.8 


9x9 




This work 


CO J = 6 - 5 


691.473076 


50 


CSO/690GHZ Rx 


10 


2.7 X 2.7 


0.6 


This work 


CO J = 7 - 6 


806.651776 


250 


HHT'^ /800GHz Rx 


10 




0.9 


Kawamura et al. (1999) 



''Half power beam width. 
''Field of view. 

•^The Heinrich Hertz Telescope. 

'^The largest angular scale detected with the interferometric observations = 16' 
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Table 2. Peak Flux Densities and Total Flux 



Wavelengths 


Telescope/ 


Shpbw^ 


Flux density'' 


Total flux'^ 


Effective radius 


Reference 


(fim) 


Instrument 


(") 


(Jy Bcam^^) 


(Jy) 


(") 




350 


CSO/SHARC II 


8.7 


59.1 (4.1) 


477 (13) 


65.5 


This work 


350 


CSO/SHARC 


11 


34 


130 


25 


Hunter et al. (2000) 


450 


CSO/SHARC II 


9.8 


27.4 (1.2) 


137 (14) 


50.4 


This work 


450 


JCMT/SCUBA 


8 


29 (0.9) 


91.1(3.9) 




Williams et al. ('2004) 


450 


JCMT/SCUBA 


12 


32.0 (6.0) 


162 (24) 




Cesaroni et al. ( 1999) 


850 


JCMT/SCUBA 


14.4 


5.6 (0.04) 


21.9 (0.5) 


15 


WiUiams et al. ("2004) 


850 


JCMT/SCUBA 


21 


5.3 (0.6) 


19 (3) 




Cesaroni et al. ( 1999) 


1200 


Pico Veleta/MAMBO 


11 


1.1 (0.3) 


5.8 


20 


Beuther et al. f2002) 


1350 


JCMT/SCUBA 


22 


1.74 (0.20) 






Cesaroni et al. ( 1999) 


2000 


JCMT/SCUBA 


33 


0.56 (0.05) 






Cesaroni et al. ( 1999) 



''Half power beam width. 

''Error is given in the parenthesis when the literatures provide. 



Table 3. IRAS 20126+4104 CO Outflow Parameters 



Transition 




1=^ 






PA 






Tdyn 


Reference 




Red 


Blue 


Red 


Blue 




Red 


Blue 








(pc) 


(pc) 


(Mq) 


(Mq) 


n 


(km s~^) 


(km s"-"^) 


(yr) 




J = 6-5 


0.23 


0.16 


13 


28 


50 


35 


30 


> 1.7 xlO"' 


This work 


J = 1 - 


0.64 


0.35 


20 


33 


-10 


25 


> 28 


6.4x10'' 


Shepherd et al. ("2000) 



^Outflow lobe radius. 

''Mass of the outflow lobes. 

'^Outer boundary of the wing components. 

"^Dynamical timescale. 



-22- 



350/xm 450/xm 350/xni + Optical 
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Fig. 1. — SHARC II 350yum (left panel) 450/im (middle panel) images of the dense clump. 
The 350/im image is overlaid on the optical image (DSS 0.5/im, color) in the right panel. 
Blue stars mark the positions of foreground stars. Grey scale contours of two images from 
left are drawn at 2a, Sex, 9cj, 27a Slcr, and 243ct, where 1(T=200 mJy beam~^. The yellow 
cross at the center represents the position of the clump center, namely the peak of the 350/xm 
continuum emission. The beams is shown in the left bottom corner of each panel. 
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Fig. 2. — Radial column density profile of the dense clump (black dots) derived using the 
350/im data. The thick broken line represents the beam shape. The vertical dotted line marks 
-Reff350- The dash-dotted and long-dashed lines mark the effective radii of the continuum 
emission measured at 1.2mm and at 850/im, respectively. 
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Fig. 3. — 350;um/450/im flux ratio map. The 450/im emission is drawn in white contours 
corresponding to 3a, 5a, 9a, 27a, and Sla. The pink dotted contours correspond to the 3a 
and 5a levels of the 350/im emission. The yellow cross marks the clump center. 
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Fig. 4. — Plot of the 350//m/450//m flux ratio as a function of the clump radius. The bars 
are the standard deviation of the ratio in each annulus. 
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Fig. 5. — (a) CO 6 — 5 and 7 — 6 spectra observed towards the center of the clump. The sohd 
and dashed bars under the spectra indicate the ranges of red- and blue-shifted components 
(see § I3.2.3p . (b) Ratio of the brightness temperatures of the CO 7 — 6 and 6 — 5 emission. 
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Fig. 6. — Comparison between CO 6 — 5 line emission close to the systemic velocity (black 
solid line contours) and the 350/im dust continuum emission (grey scale contours, same as 
Figure [1] left panel). Contours of the CO 6 — 5 map are drawn in steps of 2(T, starting from 
the 2a level (10.6 K km s"^). 
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Fig. 7. — The colored pixels in the diagram represent the LSR velocities of the CO 6 — 5 
narrow line-width components (contours, same as Figure [6]), overlaid on the map of the CO 
6 — 5 systemic velocity component (contours, same as Figure [6]). The black cross denotes 
the clump center. 
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Fig. 8. — CO 6 — 5 spectra observed towards the positions at offsets (+45", 0) (east) and 
(+100", 0) (west) from the clump center. The dotted vertical hne marks Vco6-5- 
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Fig. 9. — The grey scale represents the non-thermal velocity dispersion (cr^„J of the nar- 
row line-width components, overlaid on the map of CO 6 — 5 systemic velocity component 
(contours, same as Figure E]). The black cross marks the clump center. 
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Fig. 10. — Plot of non-thermal velocity dispersion versus distance from the clump center. 
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Fig. 11. — CO J — 6 — 5 spectra taken towards three different positions ((RAogset"; 
DeCofrset")=(+5, —5), (—10, +5), and (0, 0)) and average spectrum of all spectra taken 
within a central area of I'x 1'20". The solid and dashed bars under each spectrum show the 
ranges of red- and blue-shifted components. The vertical dashed line is drawn at Vco6-5- 
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Fig. 12. — CO J = 6 — 5 blue-shifted (—34.0 < Vlsr km s~^< —6.0, blue dashed contours) 
and red-shifted (0.0 < Vlsr km s~^< +30, red solid contours) emission superposed on the 
350yum image (same as Figure 1). Contours of the CO J = 6 — 5 components are drawn in 
steps of 3(7, starting from the 3a level (50 K km s^^). The yellow cross marks the clump 
cente r. The shocked kn ots (green dots) marked and numbered in the diagram are taken 
from IShepherd et al.l (120001 ). 
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Fig. 13. — CO 6 — 5 blue- and red-shifted (left, same as Figure [T2|) . and CO 1 — blue- 
and red-shifted (right, same as Figure 1 of Shepherd et al. 2000) maps, superposed on 
the 350/450/im flux ratio map. White dashed and red solid contours denote the blue- and 
red-shifted lobes respectively. 



